library(xlsx)

Data = read.xlsx('Proliferation 24 hr.xlsx', 1, header=TRUE, startRow=1)

set.seed(1234)
str(Data)

res.aov2 <- aov(Proliferation ~ Aggregate + Region , data = Data)
summary(res.aov2)

library(DescTools)
PostHocTest(aov(Proliferation ~ Aggregate + Region , data = Data),method = c("hsd"),conf.level = 0.90)
